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We present the observation of the decay B° — > XcoK* as well as evidence of B + — > XcoK* + , with 
an 8.9 and a 3.6 standard deviation significance, respectively, using a data sample of 454 million 
T(4S) — > BB decays collected with the BAR4R detector at the PEP-II B meson factory located at 
the Stanford Linear Accelerator Center (SLAC). The measured branching fractions are: B(B a — ► 
XcoK* ) = (1.7 ± 0.3 ± 0.2) x 10~ 4 and B(B+ -> XcoK*+) = (1.4 ± 0.5 ± 0.2) x 10~ 4 , where 
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the first quoted errors are statistical and the second are systematic. We obtain a branching fraction 
upper limit of B(B + -> XcoK* + ) < 2.1 x 10" 4 at the 90% confidence level. 



PACS numbers: 13.25.Hw 

Theoretical predictions of branching fractions and 
rate asymmetries in non-leptonic heavy-flavor meson 
decays are difficult due to our limited understand- 
ing of the process of quark hadronization. fn the 
simplest approximation, weak decays such as B — > 
J/0 K arise from the quark-level process b — > ccs through 
a current-current interaction that can be written as 
[c7 M (l — 75)c][s7 M (l — 7s)&], where 7 M are Dirac matri- 
ces ( fx — 0,1,2,3), 75 = 70717273 and c,c,s,b are quark 
spinor fields. The colorless current £7^(1 — 75)0, which 
can create the Jji\> , can also create the P-wave state X d ■ 
It cannot, however, create the XcO, Xc2 or h c states, so 
their appearance would have to be explained by a more 
complex hypothesis. A theoretical prediction can be ob- 
tained with the factorization hypothesis [H, assuming 
that the weak decay matrix element can be described as 
a product of two independent hadronic currents. Under 
the factorization hypothesis, B — > ccK^> decays are al- 
lowed when the cc pair hadronizes to J/ip, or Xfl-> 
but suppressed when the cc pair hadronizes to XcO [2]- 
In lowest-order Heavy Quark Effective Theory, the decay 
rate to the scalar \c0 is zero due to charge conjugation 
invariance 

The decay B + — ► \cqK + has been observed by Belle 
and BABAR with an average branching fraction (B) of 
(1.40±g;?jj) x 1CT 4 using X c0 decays to K+K~ or 
7r + 7r _ . This result is of the same order of magnitude 
as the branching fraction of the decay B + — > XciK + , 
(4.9±0.5) x 10 -4 [4J , and is surprisingly large given the ex- 
pectation from factorization. Using the hadronic \c0 de- 
cays, Belle has obtained an upper limit on B° — > XrfjA 
of 1.1 x 10~ 4 at 90% confidence level 0. No predictions 
are available for B decays to XcoK* , so the branching 
fraction measurement of B — > XcoK* should improve our 
understanding of the limitations of factorization and of 
models that do not rely on factorization. 

In this paper we report the first observation of B° — > 
XjpK* an d find evidence of the decay B + — > x c0 A*+ 
[6( . We identify XcO mesons through their decays to h + h~~ 
(h = K,ir), as XcO — > K + K~ and XcO - * ir + n~ have 
a higher branching fraction than the radiative decay to 
J/07 (J/0 — > l + l~ , I = fJ. or e), that was used in the pre- 
vious search for B — » XcoK* [7j. We identify A* + mesons 
through their decay to K®ir + , where A° — > 7r + 7r~, and 
A* mesons through their decay to K + ir~ . 

The data on which this analysis is based were collected 
with the BABAR detector Q at the PEP-II asymmetric- 
energy e + e~ storage ring. The BABAR detector con- 
sists of a double-sided five-layer silicon tracker, a 40-layer 
drift chamber, a Cherenkov detector, an electromagnetic 
calorimeter, and a magnet with instrumented flux return 



(IFR) consisting of layers of iron interspersed with re- 
sistive plate chambers and limited streamer tubes. The 
data sample has an integrated luminosity of 413 fb _1 
collected at the T(4S) resonance, which corresponds to 
(454 ± 5) x 10 6 BB pairs. It is assumed that the T(4S) 
decays equally to neutral and charged B meson pairs. 
In addition, 41 ftr 1 of data collected 40 MeV below the 
T(4S) resonance (off-resonance data) are used for back- 
ground studies. 

Candidate B mesons are reconstructed from five tracks 
for charged B decays and four tracks for neutral B de- 
cays, where three and four tracks, respectively, are consis- 
tent with originating from a common decay point within 
the PEP-II luminous region. Each of the tracks is re- 
quired to have a transverse momentum greater than 
50 MeV/c and an absolute momentum less than 10 GeV/c. 
The tracks are identified as either pion or kaon candi- 
dates, with protons vetoed, using Cherenkov-angle in- 
formation and ionization energy-loss rate (dS/drr) mea- 
surements. The efficiency for kaon selection is approx- 
imately 80%, including geometric acceptance, while the 
probability of misidentification of pions as kaons is below 
5% up to a laboratory momentum of 4GeV/c. Muons 
are rejected using information predominantly from the 
IFR. Furthermore, the tracks are required to fail an elec- 
tron selection based on their ratio of energy deposited 
in the calorimeter to momentum measured in the drift 
chamber, shower shape in the calorimeter, dE/dx, and 
Cherenkov-angle information. Candidate K® mesons are 
reconstructed from ir + n~ candidates, and are required to 
have a reconstructed mass within 15 MeV/c 2 of the nom- 
inal K° mass a decay vertex separated from the B + 
decay vertex with a significance of at least five standard 
deviations, a flight distance in the transverse direction 
of at least 0.3 cm and a cosine of the angle between the 
line joining the B and K® decay vertices and the K® 
momentum greater than 0.999. 

Four kinematic variables and one event-shape variable 
are used to characterize signal events. The first kinematic 
variable, AE, is the difference between the center-of-mass 
(cm.) energy of the B candidate and y/s/2, where ^/s is 
the total cm. energy. The second is the beam-energy- 
substituted mass toes = \J ( s /2 + p. ; ■ PB) 2 /Ef — p^, 
where p# is the reconstructed momentum of the B candi- 
date, and the four-momentum of its parent T(4S) in the 
laboratory frame, (E^Pi), is determined from nominal 
colliding beam parameters. The third kinematic vari- 
able is the Kir invariant mass, tok», used to identify 
K* candidates, where Kir is K®ir + or K + ir~ for K* + 
or A"* candidates, respectively. The fourth kinematic 
variable is the h + h~ invariant mass, rrihh, used to iden- 
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tify XcO candidates. Candidate B mesons are required 
to satisfy | AS | < 0.1 GeV, 5.25 < m ES < 5.29GeV/c 2 , 
0.772(0.776) < m K , < 0.992(0.996) GeV/c 2 for B+(B°) 
candidates and 3.35 < m hh < 3.50 GeV/c 2 . The event- 
shape variable is a Fisher discriminant T @ , constructed 
as a linear combination of the absolute value of the cosine 
of the angle between the B candidate momentum and the 
beam axis, the absolute value of the cosine of the angle 
between the thrust axis of the decay products of the B 
candidate and the beam axis, and the zeroth and second 
angular moments of energy flow about the thrust axis of 
the reconstructed B. 



Continuum quark production [e 1 



qq, where q 



— u,d,s,c) is the dominant source of background. It is 
suppressed using another event-shape variable, |cos0t|, 
which is the absolute value of the cosine of the angle 
9t between the thrust axis [n| of the selected B candi- 
date and the thrust axis of the rest of the event. For 
continuum background, the distribution of |cos$t| is 
strongly peaked towards 1 whereas the distribution is 
essentially flat for signal events. Therefore, the relative 
amount of continuum background is reduced by requiring 
| cos e T | < 0.9. 

Backgrounds from other B meson decays are studied 
with Monte Carlo (MC) events, using at least 10 3 times 
the number of events expected in data for specific decay 
modes that are the possible sources of background for 
this analysis. 

Potential charm contributions from B — > S(— > 
K*h~~)h + events are removed by vetoing events with a 
reconstructed K*h~ invariant mass in the range 1.83 < 
iriK'h < 1-91 GeV/c 2 . To remove background from D° 
mesons, a veto is applied to any Kit pair with an invariant 
mass in the range 1.83 < rriKn < 1.91 GeV/c 2 for each 
B — ► Xco{-^ h + h~)K* decay. Studies of MC events show 
that the largest remaining charmed backgrounds are B + 
-> U°(-> #>+7r-)7T + and B° -> D~(-+ K + tt~tt~)tt+ , 
with 12% and 10% passing the veto, respectively. Surviv- 
ing charmed events have a reconstructed D mass outside 
the veto range as a result of using a 7r or K candidate 
that is incorrectly selected from the other B decay in the 
event. 

A fraction of signal events has more than one B candi- 
date reconstructed. For those events, the candidate with 
the highest x 2 probability of the fitted B decay vertex is 
selected. Studies of MC events show that less than 11% 
of events are reconstructed from the wrong candidate, 
where these incorrectly reconstructed events are modeled 
in the fit to data. 

After applying all selection criteria, there are five main 
categories of background from B decays: two- and three- 
body decays proceeding via a D meson; non-resonant 
B — > K*h + h~ and B — > KttXcO', combinatorial back- 
ground from three unrelated particles (K*h + h~); two- or 
four-body B decays with an extra or missing particle and 
three-body decays with one or more particles misidenti- 



fied. Along with selection efficiencies obtained from MC 
simulation, existing branching fractions for these modes 
3, EH are used to estimate their background contribu- 
tions that are included separately and fixed in fits to 
data. For the non-resonant backgrounds, where there is 
no branching fraction information, fits to sideband data 
(0.996 < m K , < 1.53 GeV/c 2 and 3.2 < m hh < 3.35 
GeV/c 2 ) are performed to estimate the background con- 
tributions. 

In order to extract the signal event yield for the channel 
under study, an unbinned extended maximum likelihood 
fit is used. The likelihood function for N events is 



C 



— ex P -]>> n hr 

\ i=i / j=i \i=i 



rii Pi(d, Xj) 



(1) 



where M = 3 is the number of hypotheses (signal, contin- 
uum background, and B background), m is the number of 
events for each hypothesis determined by maximizing the 
likelihood function, and Pi(a, x~j) is a probability density 
function (PDF) with the parameters a and variables x 
= ( m ES> AS, T, rriK* and mhh)- The PDF is a product 
Pi(a,S) = Pi(d mES ,m ES ) x Pi(a A E, AS) x P^dr,^) x 
Pi{a mK , ,m K *) x P l {d mhh ,m, h h)- Studies of MC simu- 
lation show that correlations between these variables are 
small for the signal and continuum background hypothe- 
ses. However, for B background, correlations of a few 
percent are observed between toes and AS, which are 
taken into account by forming a 2-dimensional PDF for 
these variables. 

The parameters for signal and B background PDFs are 
determined from MC simulation. All continuum back- 
ground parameters are allowed to vary in the fit, in or- 
der to help reduce systematic effects from this dominant 
event type. Sideband data, defined to be in the region 
0.1 < AS < 0.3 GeV and 5.25 < to es < 5.29GeV/c 2 , 
as well as off-resonance data, are used to model the con- 
tinuum background PDFs. For the Toes PDFs, a Gaus- 
sian distribution is used for signal and a threshold func- 
tion [13] for continuum background. For the AS PDFs, 
a sum of two Gaussian distributions with distinct means 
and widths is used for the signal and a first-order polyno- 
mial for the continuum background. A two-dimensional 
(m E s, AS ) histogram is used for B background. The sig- 
nal, continuum and B background T PDFs are described 
using a sum of two Gaussian distributions with distinct 
means and widths. For ttik* PDFs, a sum of a relativis- 
tic Breit-Wigner function Q and a first-order polyno- 
mial describes each of the signal, continuum, and B back- 
ground distributions. Within the ttik* fit range, there is 
also the possibility of B background contributions from 
non-resonant and higher K* resonances; these contribu- 
tions are modeled in the fit using the LASS parameteri- 
lH , l3 | . The contribution from this background is 



zation 



estimated by extrapolating a Kir invariant mass projec- 
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tion fitted in a higher-mass region (0.996 < to_r-* < 1.53 
GeV/c 2 ) into the signal region. This estimated back- 
ground is modeled in the final fit to the signal region and 
assumes there are no interference effects between the Kn 
background and the if* (892) signal. Finally, for rrihh 
PDFs, a sum of a relativistic Breit-Wigner function and 
a first-order polynomial is used to describe signal and a 
first-order polynomial to describe the continuum and B 
background distributions. The non-resonant h + h~ back- 
ground is modeled by a first-order polynomial, and the 
background is estimated by extrapolating the invariant 
mass projection fitted in the lower mass region (3.2 < 
rrihh < 3.35 GeV/c 2 ) into the signal region. The signal 
first-order polynomial component of the ttlk* and rrihh 
PDFs is used to model misreconstructed events; for ex- 
ample where a K + from the K*° is reconstructed as a 
Xco daughter particle, and vice versa. 

To extract the B — + XcoK* branching fractions, B, the 
following equation is used: 



B 



N B b x e x B{xc0 -> h+h-) 



(2) 



where n s i g is the number of signal events fitted, e is 
the signal efficiency obtained from MC and Nbb is the 
total number of BB events. The efficiencies take into 
account both B(K*+ -> K°tt+) = 2/3 and B(K*° -> 
K + tt~) = 2/3, assuming isospin symmetry, as well as 
B(K° -> K°) = 1/2 and B(K° S -> 7r+7r-)[|. The branch- 
ing fractions are calculated taking into account B(xcO ~ ¥ 
K+K~) = (5.5 ± 0.6) x 10~ 3 and B{ X m -> tt+tt-) = 
(7.3 ±0.6) x 10" 3 0. 

We observe the decay B° — > XcqK* with an 8.9 stan- 
dard deviation significance and measure the branching 
fraction B{B° -> XcoK* ) = (1.7 ± 0.3 ± 0.2) x 10~ 4 . 
We find evidence for B + — > XcoK* + with a 3.6 stan- 
dard deviation significance and set a 90% confidence 
level upper limit on the branching fraction of 2.1 x 
10~ 4 . Figure [1] shows the fitted mgs and rrihh projec- 
tions for the B+ -> XcoK*+(xcO -> K+R-), B+ -> 
Xc K*+( Xc0 -> 7T+7T-), B° ^XcoK*°( X cO - K+K-) 
and B° — > XcoK*°(XcO — ¥ 7r + 7r~) candidates, while the 
fitted signal yields, measured branching fractions and 
upper limits are shown in Table |TJ The candidates in 
Fig. [1] are signal-enhanced, with a requirement on the 
probability ratio V s i s /(V a i g + Pbkg)j optimized to en- 
hance the visibility of potential signal, where V S ig and 
Vhkg are the signal and the total background probabil- 
ities, respectively (computed without using the variable 
plotted). Figure [5] shows the — 21n£ distributions for 
both B + — > Xc()K* + and B° — ► XcoK* as a function 
of branching fraction. The — 21n£ distributions for the 
final states (xco - * K + K~ and Xco ^ tt + tt~) are com- 
bined to give final branching fractions shown in Table [J 
The 90% confidence level branching fraction upper limit 
(Btjl) is determined by integrating the likelihood dis- 



tribution (with systematic uncertainties included) as a 
function of the branching fraction from to i3uL, so that 
J o euL CAB — 0.9 J °° CAB. The signal significance S, in 
units of standard deviation, is defined as v2A ln£, where 
A In £ represents the change in log-likelihood (with sys- 
tematic uncertainties included) between the maximum 
value and the value when the signal yield is set to zero. 
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FIG. 1: Maximum likelihood fit projections of niES (left col- 
umn) and mhh (right column) for signal-enhanced samples 
of B — » XcoK* candidates. The dashed line is the fitted 
background PDF while the solid line is the sum of the signal 
and background PDFs. The points indicate the data. The 
plot shows projections for B + — » XcqK* {XcO — * K + K~) 
(a) and (b), for B + -> XcoK* + (XcQ -> tt + tt~) (c) and (d), 
for B° -> XcoK*°( X cO -» K+ K~) (e) and (f), and B° ~* 
XcoK*°(xc0^k+tt-) (g)and (f). 



Contributions to the branching fraction systematic un- 
certainty are shown in Table [TTJ The presence of a non- 
resonant K + K~ and 7r + 7r~ can give rise to interference 
effects, resulting in a departure from the rrihh PDF used 
in the fit to data. In order to estimate how much this 
can affect the extracted yields, the fit is repeated with 
the inclusion of a PDF describing the interference be- 
tween the Breit-Wigner and non-resonant amplitudes in 
the rrihh distribution. This shape consists of the squared 
modulus of the sum of a Breit-Wigner and a constant 
amplitude, carrying an arbitrary phase difference. The 
relative weight of these two components and their phase 
difference are allowed to vary to obtain the best fit. The 
signal yields derived from this fit are larger than the nom- 
inal fit in Table Q] and the difference from the nominal 
fit is used as an estimate of the systematic error in Ta- 
ble [IT] due to neglecting interference effects. Interference 
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TABLE I: Total number of events in the fit, B background yields (B bkg), signal yields, efficiencies, and branching fractions 
£>, measured using B — > \ c oK* events. Fit bias corrections are applied to the signal yields and branching fractions. The first 
error is statistical and the second error is systematic. The significance S is shown for B + —* XcoK* + and B° — > XcoK* and 
the branching fraction upper limit, £>ul, at the 90% confidence level is shown for B + — > \ c oK* + . 



Mode 


Total 


B bkg 


Signal 


Signal 


B 


Bul S (a) 




Events 




Yield 


Efficiency(%) 


(x 10- 4 ) 


(x 10- 4 ) 


B+ X coK* + 


X co -» K+K- 


156 


8 


13 ± 5 


3.2 


1.6 ± 0.7 ± 0.2 




XcO -> + 


1065 


65 


15 ± 9 


3.8 


1.2 ± 0.7 ± 0.2 




Combined 










1.4 ± 0.5 ± 0.2 


2.1 3.6 


nU TV'* 

B -> x^oif 


Xc0 -» if + if" 


690 


20 


47 ± 10 


11.1 


1.7 ± 0.4 ± 0.2 




XcO -> 7T + 7T 


4507 


154 


72 ± 15 


12.8 


1.7 ± 0.4 ± 0.2 




Combined 










1.7 ± 0.3 ± 0.2 


8.9 



5, 

CM 
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FIG. 2: Distribution of — 21n£ as a function of branching 
fraction for B+ -> XcoK*+ (a) and B° -> XcoK* (b). In each 
case, the upper dashed line is the decay XcO — > K + K~ and 
the lower dashed line is the decay XcO — * The solid 

line is the combination of the two. In all cases systematics 
contributions are included and the — 21n£ distributions have 
been shifted vertically so the minimum value is 0. 



effects between the if* (892) and spin-0 final states (non- 
resonant and if,* (1430)) integrate to zero if the accep- 
tance of the detector and analysis is uniform; the same is 
true of the interference between the if* (892) and spin-2 
final states (if|(1430)). Studies of MC events show the 
efficiency variations are small enough to consider these 
interference effects insignificant. The integrated inter- 
ference between if* (892) and other spin-1 amplitudes 
such as if* (1410) is in principle non-zero, but in prac- 
tice is negligible due to the small branching fraction of 



TABLE II: Summary of systematic uncertainty contributions 
to the branching fraction measurements B — > XcoK* . Multi- 
plicative and additive errors are shown as a percentage of the 
branching fraction. The final row shows the total systematic 
error on the branching fractions. 



Error 




XcoK* 


f XcoK*° 


XcoK* 


Source 


Xco[KK) 


Xo0(7T7r 


) Xco(KK) 


XcoOtt) 


Multiplicative errors (%) 








Interference 


7.2 


8.3 


6.8 


10.1 


Tracking 


4.0 


4.0 


3.2 


3.2 


K° s Efficiency 


1.7 


1.7 






Particle ID 


1.9 


2.7 


2.4 


3.2 


B(xco h+h- 


10.9 


8.2 


10.9 


8.2 


No. of BB 


1.1 


1.1 


1.1 


1.1 


Tot. mult.(%) 


13.9 


12.8 


13.8 


13.8 


Additive errors 


(%) 








Fit Bias 


1.3 


4.4 


1.8 


3.9 


B background 


0.5 


4.5 


1.4 


1.5 


PDF params. 


0.6 


3.4 


0.3 


2.6 


Tot. add. (%) 


1.5 


7.2 


2.3 


4.9 


Total (10~ 4 ) 


0.2 


0.2 


0.2 


0.2 



if*(1410) -> K+ir- (6.6 ± 1.3% 4]) and the fact that 
the if7r mass lineshapes have little overlap. Errors due 
to tracking efficiency, ifg reconstruction efficiency and 
particle identification are assigned by comparing control 
channels in MC simulation and data. The branching frac- 
tion error of XcO ~ * h + h~ is taken from the combination 
of previous measurements [3| ■ The number of BB events 
is determined with an uncertainty of 1.1%. To estimate 
errors due to the fit procedure, 500 MC samples contain- 
ing the numbers of signal and continuum events mea- 
sured in data and the estimated number of exclusive B 
background events are used. The differences between the 
generated and fitted values are used to estimate small fit 
biases (see Table [TTJ) that are a consequence of correla- 
tions between fit variables. These biases are applied as 
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corrections to obtain the final signal yields, and half of 
the correction is added as a systematic uncertainty. The 
uncertainty of the B background contribution to the fit 
is estimated by varying the known branching fractions 
within their errors. Each background is varied individ- 
ually and the effect on the fitted signal yield is added 
in quadrature as a contribution to the uncertainty. The 
uncertainty due to PDF modeling is estimated by vary- 
ing the PDFs by the parameter errors. In order to take 
correlations between parameters into account, the full 
correlation matrix is used when varying the parameters. 
All PDF parameters that are originally fixed in the fit are 
then varied in turn, and each difference from the nominal 
fit is combined in quadrature and taken as a systematic 
contribution. 

In summary, we have observed the decay B — ► XcoK* 
with an 8.9 standard deviation significance and find ev- 
idence for B + — ► XcoK* + with a 3.6 standard deviation 
significance, placing an upper limit on the branching frac- 
tion. The B Q — ► XcqK* branching fraction does not 
agree with the zero value expected from the color-singlet 
current-current contribution alone, and is approximately 
half the B a -> XaK* branching fraction ((3.2 ± 0.6) x 
1CU 4 [4]), which is surprising when taking into account 
factorization expectations. 
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